Violation of a temporal Bell inequality for single spins in solid 
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Quantum non-locality has been experimentally investigated by testing different forms of Bell's 
inequality, yet a loophole-free realization has not been achieved up to now. Much less explored 
are temporal Bell inequalities, which are not subject to the locality assumption, but impose a 
constrain on the system's time-correlations. In this paper, we report on the experimental violation 
of a temporal Bell's inequality using a nitrogen vacancy defect (NV) in diamond and provide a 
novel quantitative test of quantum coherence. We therefore present a new technique to initialize 
the electronic state of the NV with high fidelity, a necessary requirement for reliable quantum 
information processing and/or the implementation of protocols for quantum metrology. 



Ever since the birth of quantum mechanics, its non- 
determinism and non-local structure have motivated in- 
creasing theoretical scrutiny. More recently, discussions 
have moved to the experimental realm but many as- 
pects of the initial debate still remain to be fully un- 
derstood. The superposition principle underpins quan- 
tum non-determinism and its application to systems of 
increasing complexity can lead to contradictions with 
our direct experience of the world, as neatly exempli- 
fied by the Schrodinger's Cat Gedankenexperiment [1]. 
Schrodinger envisaged a (quantum) cat trapped inside a 
box equipped with some device which may, randomly, kill 
the cat. Quantum mechanics tells us that at any time, if 
unobserved, the cat is both dead and alive. However, if 
the system is probed, for instance by opening the box, we 
will find the cat in one of the two possible states, either 
alive or dead. According to the formalism of quantum 
mechanics, this result is non-deterministic, i.e. the out- 
come of our measurement, and therefore the (classical) 
state of the cat is not even defined before opening the 
box. But how can we know that quantum mechanics is 
not simply incomplete, and that there are not other, "hid- 
den" variables, on which the state of the cat depends on, 
thus rendering the measurement deterministic? When 
formulated in the context of spatially separated subsys- 
tems, this type of argument leads to the celebrated Bell's 
inequality, which limits the strength of the statistical cor- 
relations between distant subsystems within a local real- 
istic theory [2H5]- 

On the other hand, in contexts where locality is not 
relevant, the focus is on the premises of realism and the 
characterization of the type of temporal (vs spatial) cor- 
relations that would emerge within that realist descrip- 
tion. This issue was first investigated by Leggett and 
Garg in 1985 ^ and led to the formulation of the so called 
temporal Bell inequalities (TBI). The major difference to 
the original Bell inequalities is that instead of correlations 
between states of two spatially separated systems, one is 



now concerned with correlations of the state of a single 
system at different points in time. In [5], a measurable 
inequality is obtained from the assumptions of macro- 
scopic realism (a macroscopic system is always in one of 
its macroscopically distinct states) and the possibility to 
perform non-invasive measurements (measurements do 
not infiuence the dynamics of the system) . It is clear that 
these assumptions, and therefore the resulting inequality, 
are typically violated by quantum mechanics. However, 
we would expect these assumptions to become valid at 
some point as the considered system becomes closer and 
closer to a truly macroscopic object. Within this view, 
TBI provide a criterion to characterize the boundary be- 
tween the quantum and classical domains and the pos- 
sible identification of macroscopic quantum coherence. 
Leggett and Garg suggested an experimental scenario for 
testing their inequality using superconducting interfero- 
metric devices (SQUIDs), whose macroscopically distinct 
states would correspond to counter-propagating currents. 
Other proposals for experimental tests followed ^i8j, but 
the actual experimental implementation was hindered by 
the assumption of non-invasive measurability, which is 
violated by projective measurements. There were nev- 
ertheless some suggestions to circumvent this by using a 
measurement strategy, which, together with the realism 
assumption, would result in non-invasive measurements, 
e.g. a "negative result measurement" as suggested by 
Leggett and Garg, or delayed-choice measurements [H]. 
A more recent approach is to use weak measurements, 
assuming that the dynamics of the system (including 
possible hidden variables) is thereby disturbed only very 
slightly. Using this method, the first experimental vio- 
lation of temporal Bell inequalities was recently demon- 
strated in superconducting qubits [3]. 

Here we will follow a different approach and use the 
temporal Bell inequality presented in [101 |TT] based upon 
the strategy of replacing the condition of non-invasive 
measurements by an additional assumption of stationar- 
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ity of the correlations. The advantage of this formula- 
tion is that it leads to easily testable inequalities (using 
projective von Neumann measurements), whose violation 
can be interpreted in a transparent way. Specifically, we 
will test a temporal Bell inequality based upon the fol- 
lowing assumptions: 

Al) "Reality": The state of any physical system is 
always well defined, i.e. the dichotomic variable 
Mi{t), which tells us whether {Mi{t) = 1) or not 
{Mi(t) = 0) the system is in state i, is, at any time, 
M,(t) = 0,1. 

A2) "Stationarity" : The conditional probability 
Qij{ti,t2) to find a system in state j at time t2, 
if it was in state i at time ti only depends on the 
time difference t2 — ii. 

Effectively, the stationarity assumption implies restric- 
tions on possible deterministic hidden variable theories, 
which cannot be ruled out by experimental violation of 
the inequality (cf. [TT] for details). However, the as- 
sumed stationarity of any two-time correlation function 
is in principle amenable to separate experimental test- 
ing and expected to hold not only for idealized closed 
quantum dynamics but also in open systems subject to 
purely Markovian noise at a rate 7, where idealized si- 
nusoidal two-time correlations would be exponentially 
damped with a factor 'j{ti — tj) [12]. In a real experi- 
ment, a non-stationary dynamics could arise as a result 
of possible non-Markovian interactions. In our case, the 
environmental back-action of the environment onto the 
system is negligibly small in the observed time scales and 
the considered parameter regimes, which endorses the va- 
lidity of the stationarity assumption in this scenario and 
therefore enhances the physical character of the hidden 
variables theories put to test. 

Following [TT], the reality and the stationarity assump- 
tions yield the constrain 

Q,,(0,2i)-Q|(0,t)>0. (1) 

This inequality bounds the strength of the temporal cor- 
relations that can arise in a classical framework and, like 
the original TBI, can be violated by quantum mechanical 
unitary dynamics, e.g. by a Rabi oscillation or Larmor 
precession. When the system's evolution is not closed, 
violation of the classical bound can persist depending 
on the noise strength. For the case of Markovian noise, 
above a certain value of the noise intensity and after cer- 
tain time, the dynamical evolution will not longer violate 
the classical limit and the observed conditional probabili- 
ties could be accounted for in terms of statistical mixtures 
of orthogonal states. On the other hand, when the classi- 
cal limit is violated, the system has a degree of coherence 
leading to observable probabilities that cannot be sim- 
ulated without resorting to superposition of orthogonal 



states, which are alien to a realist description of the dy- 
namics. Within this view point, an experimental test of 
the inequality ([T]) provides supplementary, quantitative 
information to traditional schemes to test for quantum 
coherence in terms of Ramsey fringes, given that it sets 
a maximal loss of visibility of the fringes pattern in a 
given time interval in order to violate the classical bound 
and therefore ensure the impossibility of an alternative 
hidden variables description. 

We now describe the experimental set up and the ob- 
tained results. In order to measure the conditional prob- 
abilities in the TBI ([T]), we need projective quantum non 
demolition (QND) measurements. We will use the ni- 
trogen nuclear spin associated with the nitrogen- vacancy 
defect (NV) in diamond to experimentally test this in- 
equality, using repetitive readout |13' and high magnetic 
fields that allows "single-shot" measurements [l^. An- 
other approach to QND measurements of the NV is dis- 
persive coupling to light [TS]. The basic level scheme of 
the negatively charged NV (NV~) is depicted in fig. [l]a). 
Green light brings NV^ into its excited state (ES) while 
conserving the spin polarization, from where it has two 
different decay paths depending on the electronic spin 
state. For the TO5 = state there will be predomi- 
nantly radiative decay into the ms = ground state 
(GS). On the other hand, if the spin state is ms = ±1, 
an inter-system crossing to the metastable (MS) singlet 
state is more likely to occur, from where it decays into the 
ms = ground state. Due to the (mostly) non-radiative 
nature of this decay, the fluorescence of the 7715 = ±1 
state is lower than from the 7715 = state. This dynam- 
ics allows optical spin state detection, and also polarizes 
the NV^ into the ms = ground state. 

Whereas the electron spin state is destroyed during 
this process, the nuclear spin state can be made very ro- 
bust by applying a strong magnetic field [H] {B = 0.6 
T in our case). The nuclear spin state can be mapped 
onto the electron spin with a GNOT-gate (realized by a 
nuclear spin selective tt pulse onto the electron spin, cf. 
blue arrow in fig. [l] a)), which is then optically mea- 
sured. This quantum non-demolition (QND) measure- 
ment is repeated many times (~2G00) to gather enough 
information to directly determine the nuclear spin state, 
by summing up the fluorescence photons. Fig. [l]b) shows 
a histogram of many such readouts, showing two separate 
Poissonian distributions, corresponding to nuclear spin 
state TO/ = -1-1 (left distribution) and to/ = — 1,0 (right 
distribution). By introducing a threshold between the 
two distribution, the number of collected photons can be 
directly converted into nuclear spin state. The overlap of 
the distributions introduces false results, decreasing the 
fidelity F, which is 0.91 for initialization -I- readout 

in our case. 

All these measurements are carried out on NV^. As 
has recently been shown 16], under typical measurement 
conditions the NV resides 30% of time in its neutral 
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FIG. 1: a) Energy level scheme of the NV~ defect in a small 
magnetic field. Optical transitions occur between ground 
state (GS) and excited state (ES) (vertical arrows), dotted 
lines (gray) indicate radiationless decays. Line thickness cor- 
responds to transition rates, b) Histogram of many subse- 
quent QND measurements of the nitrogen nuclear spin. Low 
fluorescence level indicates that the MW n pulse was success- 
ful, i.e. that the nuclear spin state is mi = +1. c) Magnifica- 
tion of the ms = 0, — 1 levels including hyperfine splitting due 
to the ^*A'' nuclear spin (states are denoted as \ms,mj)). The 
dotted (blue) arrow illustrates the nuclear spin selective MW 
TT pulse, the dashed (purple) arrow the nuclear spin transition 
driven by RF pulses. 



state NV*^. This is due to a two-photon-ionization pro- 
cess during illumination, at which the first photon brings 
the NV~ into its excited state, from where it is ion- 
ized to NV° by the second photon. Illumination with 
green light also restores NV~. Although QND readout 
of the nuclear spin is possible under this condition, radio- 
frequency (RF) pulses aimed for nuclear spin transitions 
in the NV" nis ~ ground state will have no effect on 
NV° (30% of times), because of different hyperfine split- 
ting. Additionally, coherent driving in NV'^ is hardly 
possible due to fast dephasing.This hinders the observa- 
tion of high contrast Rabi oscillation, necessary to show 
violation of the TBI 

In the following, we present a new technique to non- 
destructively detect and thereby initialize the charge 
state of the NV. Since the zero phonon line (ZPL) of 
NV" is 575nm, illumination of light with wavelength 
A 3> 575nm will not induce any fluorescence in this 
state. On the other hand, fluorescence of the NV~ can 
be observed for illumination wavelengths of up to 637nm, 
which is the ZPL of NV~ . The ionization rate of NV^ 
decreases quadratically with the illumination power (two 
photon process) [TS], whereas the fluorescence only de- 
creases linearly. Using low power excitation light, it is 
therefore possible to measure fluorescence from NV~ be- 
fore it becomes ionized. This allows to initialize the NV 
in its negative charge state, after the initialization of the 
nuclear spin. Fig. [2] a) shows a timetrace of the fluo- 
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FIG. 2: a) Timetrace of the fluorescence of the NV under 
cw illumination with 600nm, 0.4/iW laser light (orange line). 
The red line shows the most likely fluorescence level. The 
lifetime of NV~ during orange illumination is around 600 ms. 
b) Histogram of measurement results showing the distribution 
of counted photons during orange illumination for the inset 
sequence: Green pulses "reset" the NV and 8 ms orange pulses 
measure the charge state. 



rescence of the NV under illumination with low power 
( 0.4/iW, cw) orange Hght (A = 600nm). Fig. [2]b) shows 
a histogram of charge state measurements after a green 
laser pulse (typically used to measure the spin state of 
NV~). The two different fluorescence state correspond- 
ing to NV" (low fluorescence) and NV~ (high fluores- 
cence) are clearly distinguishable by setting a threshold 
in the middle between the two peaks. Since we only want 
to assure that the NV is in its negative charge state, this 
threshold can be shifted towards higher count rates to 
increase the fidelity. Only if we find the NV to be in its 
negative charge state, we use the result of the subsequent 
experiment. 

A nuclear Rabi oscillation obtained with this method 
is shown in fig. [s] b) . The curve is slightly shifted up- 
wards, i.e. in some times the NV is not in the ms — 
NV~ ground state (cf. [H]). This can be attributed 
to the limited fidelity of the charge state measurement, 
as well as to imperfect polarization of the electron spin, 
partly due to Ti decay during the charge state measure- 
ment. The error bars show empirical errors calculated 
from many sub ensembles of the measurements for the 
respective point. 

The measurement procedure for the Rabi oscillation 
using single-shot readout is the following (cf. fig. |3]a)): 
i) Initializing the nuclear spin by measuring its state, ii) 
Measuring the charge state of the NV, iii) The applica- 
tion of radio-frequency (RF) pulses for time r to rotate 
the nuclear spin state on the Bloch sphere, iv) Measur- 
ing the final state of the nuclear spin. This measurement 
will always yield either 1 or 0, no events are missed as 
is the case with single photon detection [T7]. Averaging 
over many such measurements directly yields the condi- 
tional probabilities (3ii(0, i), Qii(0, 2t), readily derived 
from the Rabi oscillations. The result of the calculation 
of the temporal Bell inequality [l] is shown in fig. [Sjc), 
demonstrating a clear violation of the inequality. The 
red dot in fig. [sjc), corresponding to the two red dots 
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FIG. 3: a) Measurement sequence (cf. text), b) Rabi oscilla- 
tion of the nuclear spin, introduced by RF-pulses of varying 
length T (time-axis), the blue line shows a cos-fit. Red dots 
indicate points measured with higher accuracy, c) Result of 
the temporal Bell inequality The green line is obtained 
by calculating the inequality for the fit curve of b) , the dot- 
ted black line is for comparison and shows the '"ideal"' case 
(error-free measurement of Rabi oscillation). The red dot cor- 
responds to the red dots in b). 



in the Rabi oscillation (fig. [Sjb)), gives the largest viola- 
tion of eq. [T] and was measured with very high accuracy. 
Its value is —0.209 ± 0.0039, which is over 50 standard 
deviations below 0. 

This result demonstrates for the first time, that the 
dynamics of a single nuclear spin cannot be described 
by a realist theory supplemented with hidden variables. 
In future experiments, by entangling two nuclear spin 
suitable for single-shot readout, also the contextuality of 
quantum mechanics could be experimentally tested |18j . 
For future application such as quantum information pro- 
cessing and Heisenberg limited phase estimation [111 HO] , 
refining experiments ins terms of stabilizing the NV in 
its negative charge state will be of pivotal importance. 
Otherwise, the NV resides 30% of time in its neutral 
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